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ABSTRACT 


Conjugate basis functions which are continuous through- 
out a model are used to approximate a stress field. This 
conjugate approximation has less mean error than approxima- 
tions calculated using conventional finite element methods 
and is a better approximation at extreme values of stress. 
This thesis consists of two major parts. The appropriate 
analytic expressions for conjugate stress fields are set 
forth and applied to the constant strain triangle finite 
element. Analyses are then performed with these new stress 
calculations and the results compared to some previous 


Seress calculations. 
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I. INTRODUCTION 


A. BACKGROUND 

In a recent paper Oden and Brauchli [Ref. 1] state that 
one of the difficulties encountered in using the finite 
element method to approximate the stress field in a model 
has been that the stresses calculated are at best only 
averages over each element and, in general, are discontinu- 
ous between elements. Continuous stresses can be obtained 
only if strain degrees of freedom (DOF) are introduced at 
each node. This is achieved at great computational expense. 
Some of the schemes which are used to calculate the state 
of stress at a node are mentioned in Ref. l and another 
method which is based on conjugate approximations is pro- 
posed. Іп Ref. 2 Oden expands on this method explaining 
the theory of modeling a function in general and using con- 
jugate space functions for modeling in particular. 

Briefly the methods of evaluating nodal stresses re- 
ferred to in Ref. 1 include taking an average of the 
stresses of Bach element incident at the node; using a 
weighted average method developed by Wilson [Ref. 3]; using 
averages based on equating nodal forces and element stresses 
in a procedure developed by Turner, Martin and Weikel 
[Ref. 4]; or using stiffness matrices proposed by Gallagher 
[Ref. 5] which are developed relating both the stress and 


displacement fields of each element. Among all methods, the 
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proposed method, using conjugate bases, has been shown to 


give the best fit ін a least squares sense. 


БКО SCOPE AND OBJECTIVES 
In this thesis the conjugate approximation theory is set 
forth in detail and is applied to illustrative examples. 
The appropriate analytic expressions for conjugate stress 
fields are derived for the constant strain tiangle (CST) 
finite element and incorporated into a CST computer program. 
Numerous examples of in-plane problems are worked com- 
paring both the CST solution and the conjugate basis method 
ШЕРІ) solutions and the exact solution if it is available. 
Finally, some conclusions are reached concerning the 
type of in-plane problems where the CBM may be used success- 
fully and what the cost in additional computer time may be 


fom this continuous best solution. 
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II. CONJUGATE APPROXIMATION THEORY 


A. MODELING A FUNCTION 
A function F(x) over a region or domain R can be 
modeled by another function F (x) in the domain R where 


F and R ~ R. The region R is subdivided into elements 


Bes 
where each element has two or more nodes depending on the 
degree and type of the approximation function over the 
element. Within R all elements have the same number of 
nodes and each node in the global model R is designated 

sr (A = 1,6), where G is the total number of system nodes. 
The domain of the element e is r. and each mode withinm this 
local model is designated x^, (N21, No)» where No is the 
number of nodes in an element. 


For example, a one-dimensional axial bar may have the 


following three element representation (here R = R). 





Figure 1. Modeling a Bar. 


The value of F(x) at x^ is denoted p^. that is 


Е^ = в (x5, 
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Over each element the function F(x) is represented by a 


ποσα function, f (x), which is only defined on the element 


e. Its value at xn is denoted cal that ας εὖ - £ (xD). 


Hereafter, the index e denotes element e. In the figure 
N : 
below, the functions F(x), and f, for element 2 (i.e. e-2) 


2 


are shown. 





Figure 2, Global and Local Functions. 


В. INTERPOLATION FUNCTIONS 
A set of independent local interpolation functions 

EDO, (N21, М), is defined over each element such that 
Vx GO = 1 at node N and zero at all other nodes. (For ele- 
ments with linear displacement fields, such as the elements 
used in this thesis, Ux is linear between nodes.) For a 
one-dimensional 2 node element, a set of interpolation func- 
tions for representation of a linear field is, 
Y = 1-33 J> = G (1) 

1 2 
where G is a local coordinate over the closed domain [0,1] 


and origin at local node point 1. 
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Figure 3. Local Interpolation Functions 


The local function, £ (x), is approximated over the ele- 
ment by f^ (x) where 


f“ (x) = f DCN N-1, N . (7) 


Ν 
6 6 
[The convention £ (x) = E ph (x) means summation over 


repeated indices, that is 


f (x) = fps (x) + E YS (x)t. n+ te by Os 


This summation convention will be used throughout. ] 
[Note that at a node the value of £ (x) is Ed Here 
the local function τ, is modeled by ει. linear 


Interpolation functions, рі апа 05. 
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Figure 4. Approximation to Local Function 


C. INCIDENCE OPERATORS AND ASSEMBLY 


πα ΕΝ 
The function A, and A are known as incidence operators 


N 
and are used in relating the global model R to local ele- 


ments ro. They have the value of 1 or O as follows: 


1 if node A in R is coincident with 
е, node N in r. in the assembled model. (3a) 


0 otherwise 


] if node N of element r. is coincident 
ει with node A of the assembled model. (3b) 


O otherwise 
A is called a mapping function because it gives the value 
of a local node to the incident global node. Likewise, 2 


is called a decomposition function because it gives the 


global value to the corresponding local node. 


η N "m 

F = Ау f (Note this is summed over both (4) 
е е апа М.) 

NEA 

f = QA F (5) 
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The model of a function F(x) in кетпе ο be element in- 


terpolation function can now be given by 


b e 
— c e N_A 
Elz) == ES Uy Cx) Q F E=1, N. (6) 
e=1 N=1, N 
=1, G 


where N. is the number of elements. 

Equation (63 is true everywhere except at the global 
nodes with m elements coincident. Here the value of F(x) 
is m F(x) times greater than it should be. Accordingly, at 


the nodes Equation (6) becomes 
Ë e 
= _ 1 e Ν.Δ 
F(x) = i о (7) 
e=1 = 
We note that Equation (7) gives continuity of F at 


nodal points. 


A set of global interpolation functions $, (x) are now 


defined as follows: Е е 
€ N 
) уы >) 2, except at nodes 
ф, (x) = е 
Д ге 1 
е 
1 е М 
= 2, уы OX) A at nodes ( 8) 
еш 


It will be demonstrated in the examples that >, (x) is 
a linear combination of the element interpolation functions 


are independent, the 


Us. Moreover, since the local Ux 


global ФА functions are independent also. There is a unique 
Φλ for each global node. ΙΕ the element interpolation 
functions are independent, the $4 S are independent func- 


PONS. Now 


А 


F(x) = ф(х) Е = νο εἰ (9) 


(D 
II [т] 
= 
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The domain of a global interpolation function PA extends 
over all elements incident at x^. This function has the 
value 1 at x^ and zero at all other nodal points. For ex- 


ample, with the previous linear interpolation functions, 


and the three element axial domain, we have 


Figure 5, A Global Interpolation Function. 


D. FUNDAMENTAL MATRIX 
The set of G linearly independent functions {фл} form 
the basis for a G-dimensional space 9 T G is the total 
number of global nodes. The fundamental matrix of the space 
® is defined as 
CAT = σφι,» ART 1G (10) 
The notation «QA Op? is used to denote the inner pro- 


ducet οἵ ФА and Ors Here the inner product is taken as 


<P Dp = So Oh Фр aR A,T = 1,6 (11) 
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š x A T oe 
Since ^n Фф, Фл dR 0, Сір is a positive, definite 
Matrix and its inverse can always be formed. It is defined 

as 


eee (12) 
SE “ 


This will be the transformation matrix used to form а 
new set of basis functions in the G-dimensional space.  More- 
over, since In Pr dp dR = fx фр Фл ак, CAT is also a symmet- 
ТИС matrix. 

It is convenient to construct CAT from local interpola- 


tion functions as follows. Substitute Equation (8) in 


Equation (10) to obtain 
E E 
e.. f 
u E МЧ „М е f 
map 59 >%г> 7 >. 5 | ПА Әр <%ұ» Шұ? (13) 


e : : ὃ 
Here Sa E is the inner product of the local interpola- 
tion functions. It is non zero only if e=f since each func- 
tion is only defined on its own element. The inner product 


Uy p> is called the local fundamental matrix for element 
e 
CNM’ 


Equation (13) now becomes 
(14) 


E. CONJUGATE SPACES 
A conjugate basis is one which is reciprocally related 
and interchangeable as to properties with the original 


basis. A set of G linearly independent functions on Re 
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La a so-called conjugate basis set, must satisfy the re- 
eprocity- relation, 
A A А 1 А=Г 
<ф ›фг> = åp where Op = to AZT τ» 


This means that the sets tor) and T form a countable 
biorthogonal basis of the space Ф. Any sets naim and (6,7 
which satisfy Equation (15) are referred to as conjugate 
basis of 6$. 

In order to define ΠΠ uniquely, they are formed as 
follows 


$^ (x) - c^ bp (x) (16) 


А : : 
It can be shown that the ф 5 are independent functions, 


and that ! 
n 
A Cx) ES CAT ф 62 (17) 
A s AK ae NG 
Moreover <ф ‚Hr? - «C $g Op? SMIC «Opp? 
= eo Cor = J as required by Equation (15). 


Using the incident operators defined in Equation (4), 


the local conjugate basis interpolation functions are 


e 
yr = qi e (18) 


and making the substitutions from Equations (16) and (9) 


gives 
e 
Na eee l 
Ape = UN C or =) 
= ΩΝ Az ey ΩΝ τ | E (20) 
ісі DN 


Equation (20) highlights the most important feature of the 


local conjugate basis functions. The functions for each 


а 





node of each element are linear combinations of the basis 
functions of all G nodes. This gives each local node a 
basis which is continuous over the entire model. The first 


example gives a numerical illustration of this. 


F. APPLICATION OF THEORY TO FINITE ELEMENT 


E 
Recall Equation (10) where Е(х) = F^ s X UM 2i Sen 


е=1 
alternate representation of F(x) in terms of the d basis 


is obtained as follows 


К 20 K 
E 
_ K _ М „е 
= Рф "E Ya f (21) 


Here Fr is the value of F(x) with respect to the con- 


jugate basis. It is calculated by forming the inner product 


of F and фк! From Equation (21) ме have 


- Е N Bn N 
<Е,ф,> үр <Fy¢ PR? we Fx $ S047 
_ N = 
= Fy ὃν Εκ (2:2) 


Ihe local values of F are obtained by the usual decom- 


position, 


e a E - = e 
É - FA - ων <Е,%,> = <F, Py> (23) 
E e 
СЕС *^ дары N 
EEEN тле m mn (24) 


Now the original function can be approximated in the 
conjugate basis using Equation (21). This is the procedure 
used in solving problems using the conjugate basis method 


RUBM). 
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G. BEST APPROXIMATION 
Before proceeding, it is important to verify that ap- 
Proximating a function using the conjugate basis functions 
will give the best approximation in a least square sense. 
This is proven simply and clearly with a simple example. 
Тһе function Е(х) = a is to be approximated in the in- 
terval х=0 іо х=1 by a linear polynomial F(x). Here bases 


of F are functions $171 and Ф о=х 


e 


Е = aso, + а,ф, = a, (1) + a, (x) (25) 


iuc difference between the actual function and aes 


approximation is called the error, E 


E = [F = F (26) 
The function J - J(a,) is a measure of the error over 
the interval. 


m p“ dx = TOC % а, x) - ar 
С ο 


I 


J(a,) 


Y 1 
= + 3 t S (27) 


co | ts 


2 
J(ai) a, + 8182 81 


Тһе error is minimized by setting the partial deriva- 


tion of J(a.), with respect to a,» equal to zero and solv- 


ing the resulting set of equations for a,'s. 


~ 


9Ј (а. ) 
Е _ 2 
da, = 0 = 281 + а, 3 
ου (α.) 
Е _ 1 
Ja, -. Ὁ 5 а, + 329 > (28) 
Now a, = - 1 a, = 1 and F =- 1 + x”( 29) 
1 6-2 2 6 
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Using the conjugate basis method with basis functions Wa» 


ә 


F is determined as follows: 


1. 


4. 


in this 


Assume the local basis functions are Vy = = 
апа Yo = x. 


3 1 = < 2 
Form the local fundamental matrix CNM Uy » Vy š 


obtain, 
M 

сим = 611 2) (30) 
Take the inverse 

NM m 

С ES [2 i) (31) 

NR SU 

Form the conjugate basis function Yo = ФА С Pr 


example Pr vp because there is just one element. 


pr = cla +c aq,= 4 = bx 
p? = c^l 4, « c^^ à, = - 2 + 6x (32) 
e e 
Determine a, where a, = <F, ψ.» 
1 1 1 
EM 72 ва = "m 2 el 
а, = Гох фух = 753 а, = [кх $9dx = 7 
ο ο 
| " ов i 
Approximate F(x) with F = X а, V. і-1,2 
, е=1 
i S lc. 
Е = 12 (^ 6x) + z( 2+6x) 
= - £ + x | (33) 


24 





(x) F 





Figure 6. Best Approximation 


This is the same approximate function obtained using the 
previous calculus method. Therefore, the best approximation 
of a given function F(x) in the subspace @ is the function 
F(x) whose components are the inner product of F(x) and the 
conjugate basis functions. Reference 2 contains a formal 


proof that the CBM is the best fit in the least square sense. 
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III. APPLICATIONS TO DEMONSTRATIVE EXAMPLES 


In all of the problems considered in this thesis, the 
displacement must be in the domain of the structure. For 
one-dimensional problems, the model is a bar displaced along 
the x axis. The two-dimensional model is a plane with dis- 


placements in either the x or y direction or both. 


===: 


Figure 7. One- and Two-Dimensional Domains. 


Appendix A contains a brief summary of the stress- 
strain and strain-displacement relations for the in-plane 


problem. 


A. ONE-DIMENSIONAL PROBLEMS 


Jes Nonlinear Displacement and Linear Properties 


To demonstrate certain properties of conjugate ap- 


proximation functions, example 9.2 in Ref. 2 will be 
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explored iMdetazl. 
a. Statement of Problem 
Determine the stress in a nonhomogeneous,3 unit 
length bar based on approximate displacement fields. The 


stress is given by 


g(x) = k(x) 00. (34) 
where k(x) = k (l+x) (35) 


is the material modulus in units of Тс and u(x) 


is the displacement field. Assume that 
2 
Ше αἱ ο] | (36) 


where a is a small constant and x is measured in units of 
length. 
With this displacement field, the exact stress 


distribution is 


2ak 
ο 
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o(x) = - x (1+x) (37) 


The model R consists of three one-dimensional 


elements, each of unit length, that is R = R = [0,3] = 
3 
U Q. 9 
e=] 
(1) (2) (3) 
xt х? x? x4 


Figure 8. Model for Problem иа 
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The local basis functions, for alinear cs 
corresponding to a typical element, e, are the interpolation 
functions shown with Equation (1), 

e e 
ψ1(Ε} = 1 - £ and 05(5) = Б 
b. Global Interpolation Functions 
The global interpolation functions are calcu- 


lated using Equation (9), 


} 3 ÓN e 
$ (x) = Σ όλ Yy (3) А - 1,4 
е-1 N = 1,2 
NOTE: The notation (x) will be omitted unless it is needed 
mere Clarity. 
The above expression gives, 
1 1 2. 2 3 3 
Ben Desh 1,2 2185 1.3 ОЗ 
ZEE E o UE EE IE 
= 1 Y ο жо +0 +0 +0 


i11 ізі 


2 2 
I 2 2.2 
Oo = 985041 + Mody + MOY, FLOW, + NOY) + 9505 


1 2 3 
Фа= 2] + 2,99 + 2,07 + 2,0, + 2,47 + 2,0) 
= 0 + 0 +0 +0 + 0 DEM (38) 


These are shown in Figure 9. 
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0 і 4 1 Х 





1 2 3 


Figure 9. Set of Global Interpolation Functions. 
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Recall that using Equation (8) the value of ФА 
at a node where 2 or more elements are incident is still 1. 
c. Local Fundamental Matrix 
The local fundamental matrix (LFM) is formed 


using the local basis functions as discussed in Section II-D. 


оо 


е e 
a Z > 
NM Ux? UM Ыз 


1,2 
1,3 


= ош 
ІШ! 


бур = ра (Е) ФТ СЕ)аЕ = Л (1-5) (1-8) а= T 
о о 


η! 1 355] A 
] Els 1 OT e 
а = УСЕ) 91 СЕ)а = /1(&)(1-&)4&= + 


и καὶ 2 T E 
Now Em = HM 2) (40) 


Since the elements are geometrically identical, we have 


=— 
— = 


NM NM “ым” 


он 


2 3 
ο 


(1% Global Fundamental Matrix 
The global fundamental matrix (GFM) is formed 


using Equation (14), 


Z 
= 
Ф 
HRZ 
D OM 


η 4 = 
I 
оны 


Since CAT is a 4 x 4 symmetric matrix, a total 


of ten separate calculations would be needed to determine 
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all of the elements of the CEM. The calculati#oms for c 


a 
and 622 are shown below. 
il al 15 M 
л I 21114 21211 
Са = 9191 611 9194 Spot NOV ey yt 99] €292 
2.2 2.2 2.2 2 2 
Т 122 2 2451 2 οί» 
+ QII €,,* 9194 12t 98494 6,1% 9494 655 
DNE DNE 3.3 3.3 
I 1?2 3 271 3 232 3 
* 2303 C11* 2194 δη 2 1931 Ερ1 2494 522 
Е 211-12 +0+0+0 +0+0+0+0+0  +0+0+0] 
22 
( 6 | 
1.1 i m J 1 
Ет TT ie) 21211 
С; - 2:0: с11% 2:0; с.,% 2,0, с,1% 02,0; с,, 
2.2 2.2 2.2 2.2 
ο.) 12» 5 2451 2 252 2 
к 2,0: с)1% X595 c,5t* 4595 <,ү% 2,0, с,; 
a 35 3.3 3.3 
271 172 221 232 
% 8,0; сі1% $595 εφ: 9220: с,ү% 2,0; с,, 
z 210%0%0 +2 +2 404040 4040-4040] 
T; 
= 6 - (41) 
The complete GFM and its inverse are shown 
below. 
E ο ος ο. 20. 
| e O С x ν.δ» 
ΔΤ 6 De 45 o E ЕУ 
0б 01,2 E 2 56 
(42) 


In practice, the entire set of calculations 
using incidence operators is avoided and the LFM elements 
are placed in the GFM according to the following scheme. 

Each column and each row of the GFM is asso- 
ciated with one and only one global node. Likewise, each 


column and each row of the LFM is associated with just one 
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global node. In assembling the GFM, the ij element from the 
LFM associated with local nodes i,j is added to the k,l GFM 
element, where global nodes k,l correspond to local nodes 
i,j. If several elements are incident at a global node, 
then the appropriate elements from each of the connecting 
elements’ LFM will be summed in the GFM. This procedure 
avoids numerous calculations where the incident operator is 
zero. 

e. Local Conjugate Basis Functions 

The local conjugate basis functions are calcu- 


lated using Equation (19), 


e е - 1,3 
y^ (x) = e cot dn (x) N = 1,2 
e А Г 
AS Le 
The calculation for node 1 of element 1 is 
shown below. 
uu M а. "νη 
I ON η 2 J 5 i 4 
1 T 1 
Ἱ 221 2.2 1 23 1 24 
+ I» C Ф + Ф, + Q, C >. + Q, ey, 
1 1 1! 
1 З. 3 2 Hi, 3:3 ME 
t 9. € $1 + Ф, + 2, с $4 + 2 c η 
Τ᾽ 1 П 
1 51 42 1 43 1 44 
t Q. с $1 + c Ф, + Q, с $4 + m c М (43) 
e 1 
For node 1 of element 1 QA = 1 for 07 
0 for all other 
5 . AT 1 
mmerefore, substituting for C and ФА above wy becomes 


vi E 75126 UM x 7(95 + 92) ds 2 (92 + vy) 3 $2] (44) 


where the W's in brackets are the local basis functions de- 


fitrea ii Equation (i). 
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Notice that only one row ΟΕ ΕΠΕ ο еа кене 
ΕΝ 
tion was nonzero ~ that is the row where Q, - 1. Using this 


А 
fact reduces the actual calculations in a manner similar to 
the method used in forming the GFM. 
Consider the nonzero row of the above calcula- 


tion. Factoring out the common incidence operator leaves 


four terms which are the products of the [th element of the 


A row of the GFM inverse times Pr . Recall that at a node 
Pr = 1 апа the local conjugate basis functions X = x^ be- 
come equal to the A row of the GFM inverse. This greatly 


simplifies computer programming of the conjugate basis 


method. 

Returning to the example at hand. Global node 2 
is node 2 of element 1 and node 1 of element 2. Therefore, 
1 2 1 2 

2 1 2 1 
Q = 12 = 1 and Ln = о = 0 for Á = 1803, or 4. 

Showing only the nonzero terms of the calcula- 
=ron, 

2 1 zu 


N 


£ 22^ 23 24 
Y, = Y) соф, + С Q, + C $4 * C Ф 


6 1 1 2 2 3 3 
4sl-7U T 14 (05 + v0, + Шар T 245] (45) 


А similar argument for global nodes 3 and 4 


gives 


8 


V a 


p3 = $.(29] -4QUl * рү) 414005 + 92)-7021 (46) 


I 


and 


ll 


02 = 551-101 « 201 + pe) -702 + y?) + 2693] (47) 


The local conjugate basis functions are shown 


graphically in Figure 10. Notice that each function is 
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Figure 10. Local Conjugate Basis Functions. 
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continuous over the domain, x, and that the effects of nodes 
remote from the local node being considered diminish as the 
distance from that node increases. 


f. Conjugate Nodal Components and Global 
Representation 


Calculate the conjugate nodal components of the 
element stress field according to Equation (23). 
e = е 
= < > 
Ox O(x), Yy Cx) 
Since Ux are in terms of Š, б(х) wild be transformediürco 


functions of & also, 


σ(ξ) - 2ξ-: 257  £or 0« x « 1 
= 2а la a w sm< „2 
= 2E^410£412 Nor 2 ES 3 
ak 
a factor of - — has been factored out to enable compari- 


son with the results published in Ref. 2. Now 


cT - fl (28 + 28°) (1-E)ae - 2 
o 
05 = ST (2t * 26^) (D) d£ - 2 
Ὁ 
σΊ -- 02122 ie 
o 
os = 5% (252+ 6Е+ 4) (Е)аЕ = 27 
: 
σ᾽ d DE Re { 27 
о Ξ 5 (25°+10&+12)(&)4& = 22 (48) 


Now the global stress field can be modeled using Equation 


ΓΕ, 3 a N 
σ(κ) = m" Sy Ya (x) N = 1,2 
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At the global nodes we have 





so) 2352 714 19143274. 474 59 2- 
: ШЕ ТЕ ШЕ ӨШ 15 Aus qus 
2202. 
ос ү = 
DINEM A с) 
= 3.43 (49) 


Similar calculations are made for х=2 апа τοι 
Ehe results are 


m2.) 1.1087 


O (>) 235 0 
These results are plotted on Figure 12 and axe 
also shown in Table 1. Since the exact stress field was 
known, a comparison of accuracy is made. 
е: Conventional Finite Element Solution (C.F.E.) 
Here we suppose the continuous displacement 
field U(x) is approximated by the piecewise linear field, 


U(x) where 
1 го 
U(x) = 199, (x) + 89,(x) + 264(x)]. (50) 
The $'s are shown in Figure 9. This linear 


field was constructed to give U(x) = U(x) at the nodal 


points by evaluating U(x) at the nodes. The strain is now 


do, d$, 
Jot тет ο. DNS νε. 
= аф аф 
dU(x)_ _a 2 S m 
dx 9 o que t o EE Ses NE 
— = 5 2 « x «390587 
dx == = 
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The conventional solution is obtained trom 


_ dU (x) dU(x) , 
о(х) - К(х) τς. where the value of πη 


given above. 
Zs Liñear Displacement Model and Lincar Properties 
In the first example the exact stress function was 
known and was used to calculate с (х). Consider the re- 


sults if the stress field is approximated by 


4 


E dU (x) 
o(x) = k(x) dE (52) 
where A. is a constant over each element. 


This is the constant strain approximation. Since 


a constant strain triangle (CST) finite element computer 


è e . . ә 
program is used to provide the σε 8 in the two-dimensional 
ә è dU e 
problems considered later, using E 7 constant provides a 


guide to the error which may be expected. 

U(x) as given in Equation (36) is calculated at 
each node and plotted vs x below assuming linear displace- 
ment between do í The slope of each line segment is gu 


dx 


and is shown in parentheses. 
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U (x) 


2 С-1) 
8 
" (=з) 
Do.) 
5 
op 
0 1 2 3 
Figure ll. Assumed Linear Displacement Model 


e 
The conjugate nodal components, ox GO, are now 





given by 
ak 
e с o „L du e 
Oy (x) = 9 d [(1+x) ZI vn (X) dx (53) 
The results of the integration, disregarding the 
ak 
- ΙΝ factor, are 
E 221 3 _ 50 
91 7 6 ο. ME ο 
E ES NES 
2s 7 16 Оа 6 SE 


The conjugate function representation of stress is 
given by Equation (21) and the results are plotted on 


Figure 12 and included in Table 1. 
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ΝΤ, 
24 - ak, 
o: Exact 
x 
о СЕМІ 
x 
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x 
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x 
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1 2: ў 
8 
4 
0 І 
І 2 
Figure 12. Exact Stress Distribution and Approximations 


to the Exact Stress. 
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TABLE I. Summary of Results for Problems τι А 








TT AZ. 
c (x)lio (x) Error о (ж) Error c (x) Error 
BNMExact|l-CBM1 % CBM2 2 Gel ne) eG eae Ж 
0 AAA I em Mm Cd REA > сс id utend 
0, 0 -0.33 su -0.13 ГТ 1 1 s 
Hr. 4 28,481] 14.25 42 236 6.5 2,6 & 0 
Ж 12 |1127 2.50 12.06 0. 5 DoS 12 0 


ЕШ 2А 123.6 1.06 21.46 10.58 || 20 20 1:696 7 


NOTES: 1. o€x) CBM1 used known (x) 


2. 3x) CBM2 used 5. = constant 


3. O(x) C.F.E. is the conventional solution 
shown in Section III Alh 


4. σ(κ) C.Ave is the average at a node of 
the C.F.E solutions 


3. Error 4 is Exact-Approximate |400 


Exact 
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Eigure (12) shows three different approxima ions Tomiie 
exact solution. The conjugate basis method with the exact 
stress distribution gave the best solution at node 4 and 
had the least mean square error throughout. The next best 
approximation was the conjugate basis method where the 
exact stress distribution was not known. The conventional 
finite element solution is also plotted. 
3. Linear Displacement Model and Constant Properties 

This is the same problem considered in the previous 
section only this time, k(x) = k - 1. 

Equation (53) becomes 


e as dU e 
E) SS с ) πο CONSE (56) 





The nodal stresses are 


MES ег: 
η ΙΡ 
л мез 
1 IA 
ВСИ 5 
s k oo E; (57) 


The conjugate basis representation of the stress 
field is found using Equation (21) and is shown on Figure 
13. The exact solution in this case is O(x) 9 2x. The 
conventional finite element solution is also shown. The 


ок 
factor - Bom has been disregarded. 
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Figure 13. o(x) vs x for k(x) = 1. 


arnt Orme Stearn 
Consider the same three element models with k(x) = 
600 to cancel the factors 1/40 and 1/15 resulting from the 


x ч : 
10 ' This is a case 


numerical calculations and U(x) = 
where node 4 is displaced 0.3 units along the x axis. 


dU (x š : 
4000 = constant - .l In this case the exact solution, 


the conventional finite element solution, and the conjugate 
basis solution are exactly the same. This is shown in Fig- 
ure 14. This problem corresponds to a uniform bar under a 


concentrated end load P - 30. (EA) 
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x 
Exact 
C. ENS 
30 ee | CBM 
0 X 


Figure 1á, o(x) vs x for Uniform Strain. 


D. Discontinuous Displacement 


In the previous examples, the stress was in fact 
continuous throughout the model. The CBM gave the best fit 
at the extremes and the least mean-square error over all. 
However, in problems where the stress is discontinuous due 
to a discontinuous displacement approximation, the CBM will 
not be the best fit. However, we consider examples with 
Stress discontinuity to observe how the CBM tries to mini- 
mize the total error in a least squares sense in various 
models. The direct relation between discontinuous load and 
stress obtained here is true for the in-plane and bar prob- 
lems. However, other problems in ο ies may have other 
stress-load relations. For example, for bending structures 
(beams or plates), discontinuous stress arises from discon- 


tinuous moment, and not discontinuous force. 
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a. Displacement at Node 3 
For the same three-element bar given a displace- 


ment of 0.1 units at node 3, the displacement function is 


now 
= Z 
πμ ο (58) 
20 - - 
The strain for elements one and two is a con- 
Stant 0.05 and for element three is zero. The exact stress 
| ὶ | Е dU (x) 

distribution is given by о(х) = K(x) es where k(x) = 600, 


The conjugate nodal values are computed using 


Equation (23) and are shown below. 


600 


N m ppo. 
2 0 


Q? 
Il 

αἱ 
N = 


1 2 
1 1 


Q? 
m w 
li 
Q? 


3 
zu 0 (59) 


The conjugate stress field is computed using 
Equation (21) and is shown on Figure 15 along with the 


exact and conventional solutions. 
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a, (PST) 





Figure 15. ос(х) ув х for Displacement at Node 3. 


The conjugate basis method spreads the error 
throughout the model to obtain the best least squares fit, 
thereby giving a stress field in the element where there is 
no stress. Consider the area under the stress curve as a 
measure of the average error. For the exact solution A - 
30 x 2 = 60 Banik for the conjugate basis solution A = 
60.4 units. This gives a total percentage error of 0.667. 
However, at node 3, the CBM value is only 53.3% of the maxi- 
mum stress at that point. Clearly the conventional finite 
element method wae a better approximation. 

b. Displacement at Node 2 

If the model had been displaced from node 2 

instead of node 3, the results would be similar in that the 


conventional constant strain element would be a better 
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approximation and the CBM spreads the error over the entire 


model. The results are shown below. 


μα 
Χ 





Figure 16. o(x) vs x for Displacement at Node >, 


Again comparing the areas under the stress 
Елес А (е жасы): = 30: А(СВМ) - 30.89, The emror is 2.964 


overall and 53.33% less than the maximum at node 2. 


B. TWO-DIMENSIONAL IN-PLANE PROBLEMS 

The development of the computer program which was used 
to solve the two-dimensional in-plane problems is discussed 
in Section IV. Appendix B gives a ДЫ ip tion of the dif- 
ferent finite elements which were used to obtain results to 
compare with the results obtained by the CBM. The computer 


program itself is listed in Appendix C. 
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Ihe basic model used to verify the computer program was 
а 10 х 10 x 1 plane divided into either four or eight ele- 
ments as shown below. The units used for numerical calcula- 
tions are inches and pounds. The figure also shows the 


global node numbers and the element numbering system. 


2 10" 1 | 
| | 3 2 l 


10" 
14 
6 
5 
3 Model A 4 9 8 7 
Model B 
Figure 17. Typical Two-dimensional Models. 


Each node has two degrees of freedom (DOF), which are 
in-plane displacements u and v, and for a typical node n 


they are numbered as shown below. 
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DOF(2n) = νά. 0) 


DOF (2n-1)= u = ας 25 





Figure 18. Reference System. 


The following information is pertinent to all of the 


problems. 


E 


30 x 109 ІН ЕЕЕ 
н = 0.3 


All stresses plotted are o (x,y) х 104 DENE 
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1; Model A with Uniform Edge Displacement 


329 . 12232: 
| u(0,y) = 0. 
| у(0,у) - 0. 
u(l0., y) = 1. 
v(10.,y) = 0. 





329 3 


Figure 19, σ (251 х 10°) for Problem III Bi. 


The plane is prevented from displacement in the v 
direction; therefore Εν = 0. Equation (A-1) shows that the 
Stress will be uniform across the model since the strain is 
uniform. This is the result that the CST solution gives. 
Because the CST solution is continuous for all elements, the 
conjugate basis solution cannot improve upon it. Both solu- 
tions in this case are exact, for plane stress. 

This problem is analogous to the three element bar 
which was displaced at node 4. See Section III A4. In the 
following figure, the stresses obtained by the conventional 


CST finite element solution are given in each element, while 
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the stresses obtained by the conjugate basis method (CBM) 
are the values given at the nodal points. 
2. Model B with Uniform Edge Displacement 

This problem is analogous to both the problem of 
the bar pulled uniformly at its end (III A4) and the first 
Model A plane problem (III Bl). It confirms that mesh A 
was sufficient for the problem and that the CST is a 
sophisticated enough element for this type of problem. The 
difference in boundary conditions between this problem and 
problem III Bl accounts for the difference in stress values. 
In this case Б. < 0, that is contraction in the v direction 


is allowed, thereby reducing Sis 


300 300 300 


u(0.,y) - 0. 


(LOA) I; 





380 300 300 


Figure 20. Oe (PSITX 10^) for Problem III B2. 
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I Model B with Uniform Edge Displacement, Modified 


Analytical values of о were calculated for x - 
x(10,y), x = x(5,y)and x = x(0,y)using the plane stress 
ояз (A-1). The values of strain were approximated by 
dividing the total contraction in the v direction for a 
given y by the original length 10. The strain ES was con- 
sidered constant along a constant y. The displacements 
were calculated in the CST routine. 

The values of б. calculated were 


g (0.) = 329 o (5.) - 292  , o (10.) - 300 


For x = 0 and x = 5 the CBM gives the better approxi- 


mation. At x 10 the difference in percentage of error 
between the CST solution and the CBM is less than 0.67%. 
It is also interesting to note that the average of the CST 
values at nodes 2 or 8 is the same as the CBM value. This 


gives mathematical substance to the engineering practice 


of using the average value at a node. 


u(0.,y) 
v(0.,y) 
u(10.,y) = 1 





309 
Figure 21. σ. (PSI x 107) for Problem III B3. 
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4. Model A with Linear Edge Displacement 


Here the right edge is given a linear displacement 
as shown below. Using the analytical procedure discussed 
in Problem III B3, the values of с range from 291 OR 
along the edge (x,10) to zero along the edge (x,0). The 
CBM values are the best continuous approximation available 


mor this model. A plot of σ (0) vs y along (0,y) shows 


how the CBM is better in approximating the maximum value. 


> 305 
| 
| 
| 
| u(0 Уу = 
v(0- ο) ο 
| u(10.;,y) 0 
| ӘСІП) 0 
| 
| 
| 
Dl 
24 i | 24 


Figure 22. σ (151 x TOR for Problem III B4. 
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o Exact 
X 


Os o CBM 
3001 X 
G ο ὃς 
at nodes 
200. 
G CS 
Ж 
TOO. 
Y 
0 > 10 


eave iene. see o. (0) vs y for Problem III B4. 


À comparison of the percentage of error at (x,y) 
ШИ 10.) that the CBM value is 7.42 too low while the CST 
average value is 28.52 low and the minimum CST value is 
502 low. 

5. Model B with Linear Edge Displacement 

The applied displacements in this problem make it 
similar to Problem III B4; however, the other boundary con- 
ditions in this problem allow for movement in the v direc- 
tion. In part the displacement solution of the CST program 
shows that at x=0., x=5., and x=10., Ey < 0. Therefore, 


along the bottom edge where ex is very small, the effect of 


> becomes dominant and a negative > results. 


5/3 





From the results of the CST displacement solution, 


the value of omo SY along the edge (x,0.) is 0.0026 and 


the value of E - 2% - 0.016 along x-10. and 0.015 along 


(5,y) and (0,y). Using these values of E and the average 
of the εν 5 in Equation (A-1) for stress yields 


ECHO Wir Vor 


The CST solution is attempting to reflect this in 
elements 5 and 8. The CBM does show the compressive 
stresses and of course is continuous throughout the model. 
The average of the CST values at nodes where more than one 
element are er is close to the CBM value at that 
node. For example, at node 5 the CST average is 150 and 
the CBM value is 148. This again confirms the existing 


practice of averaging. 


394 238 398 
| 
( 
[ 
| u(0.,y) = 0, 
| u(10.,y) = y/10 
“py 152 
| 
| 
| 
| 
“187 66 -83 


Mifsure 24. 0 (PST x 10^) for Problem III B5. 
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6. Model B With Piecewise Linear Edge Displacement 


This problem has boundary conditions similar to 
problem III B 4; but, since it has two elements along 
x = x(10,y) and the applied displacement is only at node 1, 
the displacement profile of the edge is only piecewise 
linear with discontinuous slope. 

The expected stress profile is greatest at (x,10) 
. where E is maximum and compressive at (x,0) where e = 0 
but a compressive E exists. Likewise ES is higher along 
(0,y) because ΕΠ πι επ»; 

The CBM values conform with the expected distribu- 
tion and give a continuous stress field which is the best 


approximation in this case. 


u(0 y) - 0” 
У СОС У) πο ο, 
u(l0.:,0.)= "0; 
ντο, ο --ο. 


ч(10.,‚10.)=1. 





Figure 293: 0, (PSI x 10^ for Problem III B6. 
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7. Model A With Center Node Displacement 


For this problem the stress is discontinuous due to 
the displacement of node 5 in the interior of the model. 

The CST solution gives an exact representation of the stress 
field. Element 1 is in compression and element 3 is in 
tension. 

The CBM values are exact at the nodes. A stress 
plane constructed linearly thru the nodes would be in error 
at all locations other than the nodes. 

If a three dimensional plot is made with stress as the 
vertical axis, the volume under the CST plot would equal 


the volume under the CBM plot at 25 x —Ó 


u(0. =) 
М (Оуу) 
ТӨП уу) 
v(10. y) 
ο ο ο) 


lI 


O O O O O 


|| 





329 -329 


Figüre 26. σι (PSI x 10%) for Problem 111 B7. 
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8. Model À With Center Node Displacement Шопен 


The model is subject to a discontinuous stress field 
due to the internal displacement of the cneter node. If 
the model allowed more flexibility, one would expect the 
stress to be discontinuous in the region of node 5 and 
continuous along the edges. The CBM does not show the 
abrupt change at node 5 but does approximate the solution 


expected of a more flexible and hence more accurate model. 


329 -66 
Ф 
0...) 0. 
v(0.,y) = 0, 
u45.5.)- πι» 
329 -66 


Rua pure 2 o. (PSI x 10^) for Problem III B8 
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9. Model B With Center Node Displacement 


This problem is similar to the previous problem on 
model A with the exception that the boundary conditions on 
E problem allow displacements in the v direction. The 
stress is discontinuous in the region around the center 
Ende, node 5. However’ it is unlikely that the stress is 
as severely discontinuous as the CST solution indicates in 
this region. Intuitively om has a maximum tensile value at 
(0,5) and a maximum compressive value at (10,5) and decreases 
in absolute value away from the line of action of the 
displacement. 

For this model, subject to a point discontinuity, 
the CBM gives a better overall representation of the stress 


field. 


60 114 47 
2 
144 
u(0.,y) - 0. 
τν VO Al; 
40 174 
: | -112 
Í 
312 -24 
144 


60 114 47 


Figure 28. б. (PSI x 10%) for Problem III B9 
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10. Model B With Uniform Centerline Displacement 


This is another problem with a discontinuous stress 
resulting from an internal displacement. Note that to 
minimize deflection in the v direction, the boundary condi- 
tion v(5,5) = 0 has been imposed. 

The expected value of o (x) for x>5 is zero. "The 
CST solution is superior from two standpoints. It reflects 
the stress discontinuity very well and it gives a better 
approximation to the stress field over the right half of 
the model. It is interesting to note that the average CST 


value at nodes 2 and 8 is the CBM value. 


Eau 


UCI R: 
u (0.2 y) 
vn.) 


N 
O O 


-109 





Figure 29. o (PSI x 10°) for Problem III B10 
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ao Model B Modified, With Uniform Center mE Displacement 


This mesh is a modification of model B and was 
designed to determine how rapidly the CBM stress would con- 
verge to zero. Each node along x-5 and x-10 was displaced 
an equal amount to force rigid body motion. 

The CST model responds to the discontinuity quickly 
as expected. The CBM does try to converge and is a better 
approximation over this part of the model than the approxi- 
mation with the previous mesh (III.B.10). The conclusion 
is that if a discontinuity occurs across an entire model, 
and there is sufficient room for the continuous function to 
correct itself, and the area of interest is not at the 
discontinuity, the CBM will be the best approximation for 


a continuous field. 


ООО 
u(5.,y) 
u(0. y) 


| I 
< [= 





317 187 -22 37 


Figure 30. D (PSI x 10%) for Problem III Bil 
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2:2: Model B Tapered With Uniform Edge Displacement 


This model is another modification to model B. The 
length has been doubled and a taper introduced by making the 
right edge one half the height of the left edge. Symmetry 
about the horizontal centerline has been maintained. The 
model simulates a tapered plane and is the two dimensional 
analog of the one dimensional tapered bar. Both the CST 
solution and the CBM give the expected stress distribution. 
However the CBM is continuous over the entire model and 


therefore is a better approximation. 


IS 


u(20.3y)= A 
m , y) 
CO Y 


| 
о о 
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Figure 31. “ΒΒ (PSI x 10°) for Problem III ри 
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13. Composite Material Models With πα ο πι 
Edge Displacement 


This problem consists of using the CBM with a CST 
input on the composite material models used by Lin, Salinas, 
and Ito [Ref. 6] and comparing the results for stress 
convergence and computer time usage. The models used in 
Ref. 6 are shown in figures (32) and (33). The PLISOP model 
shown in figure 34 has been used by Salinas subsequent to 
publishing Ref. 6. The stresses resulting from using the 
CST and LST routines on figures (32) and (33) and the PLISOP 
routine in figure (34) (see Appendix B) are shown on figures 
1555156), and (37) along with the CBM values. 

In figure (37), along the edge (5.2,y), where there 
is only one material present, the CBM yields a smooth curve. 
However, where two materials are present and there is a 
sharp change in stress due to the different Young's Modulus 
of the two materials, the CBM tries to smooth out the change 
in a manner similar to the discontinuous problems shown 
earlier. As a result several CBM model values in the region 
ОГ the change in material are not in agreement with values 
obtained from the other solutions. The results are in good 
agreement away from this area of change. 

An analysis of the amount of computer time each 


program used is shown in Table II. 
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Figure УЛ" Finite Element Grid System, ΕΠΟΣ. 
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Figure 33. Finite Element Grid System 242, DOF. 
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IV. THE COMPUTER PROGRAM DEVELOPMENT 


In developing a computer program to determine stresses 
utilizing the conjugate basis method, four major requirements 
had to be met: 

i. The inverse of the global fundamental matrix is needed. 
This requires forming the element fundamental matrices 
and merging them into the global fundamental matrix 


itself. 


ii.  An'approximation to the stress field in the element is 
required. 


iii. The conjugate nodal stress values for each element are 
formed. 
iv. The conjugate basis approximation to the stress field 


must be assembled. 

There are many existing programs which satisfy the 
second requirement. In this thesis a program utilizing the 
CST written for use in Ref. 6 has been used. It was chosen 
because it was available; it uses a simple two-dimensional 
element — the three noded triangle; and problem solutions 
using this element were available for comparison as to the 
computer time used and the accuracy of results. It was not 
an objective of the thesis to develop an optimum conventional 
finite element solution. Any other CST program or a more 
sophisticated program could have been used. The advantage 
of the CBM is that it makes the stress at a node a function 
of the entire model and not just the element. 

Therefore the development of the CST program itself will 
not be discussed with the exception of a brief comment on 


the constitutive relation matrix, C. For plane stress 
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х2 yz 
1 u 0 
E = T 1 0 
l-u 
0 0 .5(1-1) 


where {0} = [с] {е} š 


For the CST, each displacement field is approximated by 
a linear field over each element, hence the local basis 
functions are area coordinates. Тһе local fundamental 
matrix, which is the inner product of each of the basis 
functions, is therefore dependent on the total area of each 
element. The integration προ ον Felippam be ο 


and the resulting matrix is 


μι CE 
А 
12 1022 1 
ic 1 


After the area of each element is calculated in the 
normal CST MEE: the element fundamental matrix — EFM — 
is formed. It is immediately merged into the global funda- 
mental matrix — GFM. The merging routine is described in 
tion III.A.l.d. 

After the entire GFM is formed, it is inverted using 
subroutine SYINV. This is the most time consuming addition 
to the CST program. See Table II. Notice that for the 92 
element nodel SYINV used 20 seconds but for the 200 element 


model the time was 133 seconds. 
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The conjugate basis nodal stresses are calculated in 
two steps at the completion of the standard CST program. 
The assumption is made that stresses will only be calculate 
at nodes. Since ФА = ] at a node (see section II.C. and 
section III.A.1.e.) the values of the local conjugate basis 
function at a node are just elements of the GFM inverse, 
in accordance with equation (44). This assumption is not a 
serious limitation because the finite element mesh can be 
arranged to have a node at any point of interest or the 
stress can be calculated by interpolating between neighboring 
nodes. However it greatly simplifies the computer 
calculations. 

Recall that the conjugate nodal values of stress for 
each element are given by Eq. (23) as 

L ΄ 


ση - ) o (x) Vs (x) dx 


The CST value of (x) is constant for each element and is 


L 
taken from the CST solution. The solution of І Way (x) dx 
0 


fPomtound@iay Ref. 7 ad ie Uu for each N. Therefore the 
conjugate value for each mode of element e is BC g (e), 
e _ A(e) 
gr ix с (е) (60) 


The conjugate function representation of stress is given 


by Eq. (21) as 


ο ο. 


II —3t=3 


е N 
On Y (x) 


e=] 
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In the computer program 0(x) is calculated one node at 
a time. The value of ση is multiplied by the element re 
GFM inverse in the row corresponding to the global node that 
the local node N is incident to and the column corresponding 
to the global node being calculated. This multiplication 
is carried out for each element and the resulting sum is 
printed out as the nodal stress. 

The modification of an existing conventional finite 
element program to perform conjugate basis calculation of 
stress requires the. addition of three simple algorithns. 

i. An algorithm which forms the global fundamental 

matrix (GFM). This algorithm may be located in the 
element stiffness DO loop, where the element funda- 
mental matrices (EFM) are individually formed and 
whose elements are inserted in the appropriate place 
in the global fundamental matrix, Care 

ii. An algorithm to invert the global fundamental matrix 
Car and obtain ca (see equation (12). CAT is a banded 
symmetric matrix and 2222. algorithms are available. 
(In this work E. Wilson's SYINV subroutine was used.) 
This algorithm is located after Cap is formed, and is 
the most time consuming aspect of the conjugate basis 
method. 

iii. An algorithm forthe calculation of nodal stresses 
associated with equations of the conjugate basis 
method. This algorithm uses the conventional stresses 


of the finite element method and must therefore follow 


the conventional stress calculation routine. 





A listing of a CST finite element routine is given in 
Appendix C with the CBM modifications shown enclosed in 


boxes. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


As a result of this study several conclusions can be 
made Meri the use of conjugate basis functions with 
in-plane problems. 

Where the model is not subject to interval discontinu- 
ities, the CBM gives the most accurate approximation to the 
stress over the entire model for a given displacement. 

If there is a point or localized discontinuity, the CBM 
will be a poor representation of the stress near the dis- 
turbance but will be a good approximation away from the 
disturbance. Mesh refinement can localize the disturbed 
area. It is understood that points of stress discontinuity 
are usually points of interest and that mesh refinements 
generally increase computer time, therefore the use of the 
CBM ís limited. 

If there is a stress discontinuity across the entire 
model, the CBM is a poor choice if the model is too small 
or too coarse to allow the CBM to settle out. Refining the 
mesh may be worthwhile if a continuous approximation to the 
stress is — in regions away from the discontinuity. 
However the additional computer time needed places a 
premium on this marginal gain. 

If a large number of problems are run on the same model 
and the global fundamental matrix can be inverted once and 
stored, the additional computer time per problem for the 


CBM solution is minimal. The advantage of a continuous 
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stress solution can be had without using the more sophisti- 


cated elements with nodal strain degrees of freedom. 


Finally, at nodes where the stress should be continuous 


and several elements were incident at that node, the average 


of the CST stress values at that node agreed very closely 


with the CBM stress value. This gives mathematical sub- 


stantiation to an engineering practice. 


It is recommended that 

The CBM be used for in plane problems without interval 
discontinuities because the resulting stress field will 
be continuous, which the actual stress field is, and the 
stresses will be consistent with the displacements. 
Further studies with the CBM should be carried out on 
other classes of problems. For example, plate problems 
or beam problems might be better approximated with CBM 
solutions. 

Additional studies be carried out to determine the cost 
effectiveness of incorporating and using the conjugate 
basis method in existing conventional finite element 
programs. 

Ат! D be undertaken to determine if the CBM 
can be modified to yield stress discontinuity when 
appropriate. 

A study be made to establish the relation of domain 


modeling to CBM analysis. 
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APPENDIX A 


THE IN-PLANE PROBLEM 


In this discussion a plane is defined as a body bounded 
by two surfaces of zero curvature whose thickness, h, is 
mush less than its surface dimensions. For in-plane 
problems only axial loads N(x,y) are considered. There are 
no lateral loads or moments. Displacement boundary condi- 
tions can be of the type аи + Ъу = с where a, b, and c 
are specified coefficients; and u and v are displacements. 


In plane forces may also be specified on portions of the 


boundary where displacement is not specified. 


P d 





Figure A-1. The In-Plane Problem 


= 


Since k is small in comparison with the plane dimensions, 


: f 
62 << s and E ons ny << E This suggests the use o 


plane stress conditions, that is, ο ον Fass are considered 


ко Бе zero. 


The basic stress strain relations become the plane stress 


equations. 
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С = (cop с.) 
x т x y 
E | 
G = — (ue +e) 
- E 
T xy 1 -- Ἡ Exy (A-1) 
The strain displacement relations are 
_ ðu 
ade 
Ир ον 
Жаа; 
ы du ду 
ΤΗ» πα, 
1 gu ду 
Exy 222 ( ду ш Әх ) с 


Substituting Eqns. (A-2) in Eqns. (A-1), the partial 


differential equations for the-in plane problem become 








E du ду 
σ = — + ws 
x l-u? L ðX » ду J 
EE du, ὃν 
e E l-u [оч уз V ду J 
E du dv 
Xy 2(1 + u) [ ду д ox 1 27 





Stress resultants N , N , N 
= y 


having the units of force 


Xy 


per unit length (f£ a are defined as 


h/2 1/2 
f = dz N, = | f E dz 
-h/2 -h/2 
(A-4) 


Ihese last expressions imply internal discontinuity in 


stress is associated with a discontinuity in applied internal 


in plane-force. 


For plane stress the equilibrium equations are 


) F,-0 


E -o 
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IN ON, N 

quem + — = 0 (А-5) 
ΟΝ ΟΝ 

a (A-6) 
y ðX 


The finite element method yields a system of algebraic 


equations which are obtained from the above partial 


differential equations by the usual method as given by 


Zienkiewicz [Ref. 


8]. 
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APPENDIX B 


IHE FINITE ELEMENT METHOD 


An extensive description of the finite element method 
can be found in Zienkiewicz [Ref. 8]. However to facilitate 
understanding of this project, certain features of the 
elements will be presented. 

In all finite element programs certain rules must be 
followed in selecting the assumed displacement field. 

l. Tne finite element displacement field must accommo- 
date rigid body displacements. 

2. At common boundaries two elements remain in 
continuous contact with each other. Neither gaps nor over- 
lappings are permitted. 

3. Constant straining patterns must be included in 


higher order descriptions of displacement. 


CONSTANT STRAIN TRIANGLE (CST) 

For the CST EC u and v in-plane displacements are 
approximated by a linear field over each element. Since 
strain is E first derivative of displacement, this gives 


a constant strain element. 
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The area coordinates used in the program are defined as 





Figure B-1. Area Coordinates 


With these coordinates the linear displacement field 


over the element is simply 


1 
Е 
с 
+ 
ES 
c 
+ 
ES 
c 


u(x,y) = 1 


I 
p 
< 
фе 
г" 
< 
of 
= 
< 


v(x,y) 5 1 (B-1) 
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LINEAR STRAIN TRIANGLE (LST) 





Figure B-2. Linear Strain Triangle 


The LST has three nodes on each side and each displace- 
ment field is taken as a 2nd order polynomial over the 
element domain. The displacements u and v for each element 


are as follows 


1 


u(x,y) = Nu + Мо, + Мо, + Мо, + Ма + Ма 


1 2 2 33 4 4 Sa 6 6 
ρα) 
у(х,у) - МУ, + NV, t М.у. + N.V, + Nev, + Neve 
where Ny = 2L, -1 for a typical corner node 
N, = 4L,L, for a typical mid side node 


The N, are called interpolation functions for the 
assumed displacement field. Differentiating the displacement 
fields gives linear strain over each element, hence the name 


linear strain triangle (LST). 
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PLISOP 

The PLISOP program was coded by Professor Cilles Cantimn 
and his students at the Naval Postgraduate School, Monterey. 
It performs a plane stress or a plane strain analysis by 
the finite element method using numerically integrated 
isoparametric elements. These elements have four sides, 
each of which may be either curved or straight. The same 
shape functions which define the geometry also are used to 
describe the displacement. The strains will be a function 
of the number of nodes per side of the element. The 
comparisons made in this thesis were to a cubic element. 
This isoparametric element uses a cubic displacement field 
and yields parabolic strain (stress) fields over each 


element. 


Figure B-3. 12 Node PLISOP Element 


All of these elements, the CST, LST and PLISOP, with 
displacements being the only degrees of freedom,yield 
discontinuous stresses along element interfaces. The CBM 


gives continuous stress along element boundaries. 
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For each of the conventional finite elements, the order 
of the strain approximation is one less than the order of 
the displacement approximation. However for the CBM the 
displacement approximation and the strain approximation are 


of the same order. 
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